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In this theoretical study, modulation techniques are developed to support the Active Sensing of CO 2 
Emissions over Nights, Days, and Seasons (ASCENDS) mission. A continuous wave (CW) lidar system 
using sine waves modulated by maximum length (ML) pseudo-noise (PN) codes is described for making 
simultaneous online/offline differential absorption measurements. Amplitude and phase-shift keying 
(PSK) modulated intensity modulation (IM) carriers, in addition to a hybrid-pulse technique are 
investigated, which exhibit optimal autocorrelation properties. A method is presented to bandwidth limit 
the ML sequence based on a filter implemented in terms of Jacobi theta functions, which does not 
significantly degrade the resolution or introduce sidelobes as a means of reducing aliasing and IM carrier 
bandwidth. 
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1. Introduction 

NASA Langley Research Center in collaboration 
with ITT Exelis has been evaluating continuous wave 
(CW) laser absorption spectrometer (LAS) systems 
as a means of retrieving atmospheric CO 2 column 
amounts from space to support the Active Sensing 
of C0 2 Emissions over Nights, Days, and Seasons 
(ASCENDS) mission. The U.S. National Research 
Council’s report, entitled Earth Science and 
Applications from Space: National Imperatives for 
the Next Decade and Beyond, identified the need for 
an active lidar space mission. This report targeted 
a need for measurements of column C0 2 to gain a bet- 
ter understanding of C0 2 sources and sinks through 
the active sensing of ASCENDS [1]. The primary 
objective of ASCENDS is to make C0 2 column-mixing 
ratio measurements during day and night over all lat- 
itudes and all seasons and in the presence of thin or 
scattered clouds. These measurements would be used 
to significantly reduce the uncertainties in global 
estimates of C0 2 sources and sinks, to provide an 


increased understanding of the connection between 
climate and atmospheric C0 2 , to improve climate 
models, and to close the carbon budget for improved 
carbon and chemical weather forecasting and as well 
as implementing policy decisions [1]. 

The critical component in the ASCENDS mission 
is the LAS for C0 2 column measurements with 
weighting toward the mid- to lower troposphere 
using the integrated path differential absorption 
(IPDA) technique from the surface and cloud back- 
scattered returns. 

IPDA lidar systems measure the difference of the 
total gas absorption along the path length between 
two or more laser wavelengths in the path of the 
laser to a scattering target. For the Earth’s climate 
system, there are large changes in the atmospheric 
variables such as cloud and aerosol distributions, 
and the profiles of water vapor and other absorbing 
gases from regional to global scales. Also, the surface 
reflectance varies with many factors and variables, 
such as surface type, roughness, slope, vegetation, 
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and soil moisture. The variations in these atmos- 
pheric and surface variables can generate large 
changes in received lidar powers from the surface 
or other targets, which can introduce large biases 
and artifacts in determination of C0 2 column 
amounts when only a single laser with a single wave- 
length is utilized. The advantage of IPDA LAS 
system is that the strong effects of the variations 
in atmospheric states and surface reflection on lidar 
return powers would be minimized from the ratio of 
the simultaneously received lidar powers from the 
two closely spaced online and offline IPDA wave- 
lengths. All the influences on the lidar return signal 
from these environmental variables are the same for 
the two or more closely spaced wavelengths. The 
ratio tends to minimize the environmental effects 
except the differential absorption optical depth at 
the desired LAS wavelengths. Thus, with all kinds 
of the considerations including experimental factors, 
atmospheric temperature and pressure profiles, and 
range estimations, high accuracy C0 2 column mea- 
surements can be achieved from the C0 2 differential 
optical depth. 

For the purposes of this paper, we consider the two- 
wavelength approach, which uses one laser with the 
wavelength tuned to the C0 2 gas absorption line 
(online), and the other wavelength tuned away from 
the absorption feature but still within 100 pm from 
that same absorption line (offline). The differential 
absorption optical depth for the C0 2 column may 
be derived from the ratio of the received optical 
power from each wavelength. There are several dif- 
ferent approaches to making this measurement. 

C0 2 IPDA LAS measurement approaches include 
pulsed [2,3] and CW systems. Direct-detection 
pulsed systems use a technique where measure- 
ments for each of the IPDA wavelengths are sampled 
successively via sequential laser shots sampled 
through different atmospheric columns and surface 
spots at the surface. The CW approach described 
here transmits the on- and offline lidar signals to 
the surface simultaneously from a single laser trans- 
mitter, thereby allowing for simultaneous sampling 
of the atmospheric column background noise and sur- 
face reflectivity within the laser spot, thus reducing 
the potential errors. In addition, a major advantage 
in this approach is the significant reduction of 
low-frequency variations in background radiation 
and colored noise, which are filtered out in a natural 
way with IM carrier-based matched-filter ap- 
proaches by effectively limiting the detection fre- 
quency and bandwidth. This is done through the 
demodulation process shown in Fig. 2, which limits 
the frequency range to a band of frequencies about 
the carrier frequency, thus filtering out unwanted 
low-frequency noise. Within the CW approach, IM 
techniques are used to separate the simultaneous re- 
turn signals from the different IPDA wavelengths, 
provide ranging capability and discriminate surface 
returns from other atmospheric returns, such as 
aerosols and clouds. The IM technique is also useful 


for extracting the laser signals from solar back- 
ground noise, which occurs as a by-product of the 
matched-filter technique used to retrieve the science 
signal of interest. Potential IM-CW techniques 
include PN code CW [4,5], linear swept sine wave 
CW [6—11], unswept sine wave CW [12—16], nonlin- 
ear swept sine wave CW [17], etc. The baseline tech- 
nique used by the NASA Langley Research Center in 
partnership with ITT Exelis for the ASCENDS 
mission is the linear swept frequency CW [11], The 
utility of this technique for both range and C0 2 
column absorption measurements has been demon- 
strated over the past several years through 
ground-based and airborne demonstrations [11], 

A drawback of various modulation techniques, 
such as the linear swept frequency, is the potential 
for generation of sidelobes in the autocorrelation 
function, which could introduce small bias errors 
when intermediate scatterers such as thin clouds 
are close to the measurement target. This is the main 
advantage of digital techniques such as PN codes and 
nonlinear frequency modulation. Previously, PN co- 
des have been described using a pure maximum 
length (ML) sequence with time shifting to separate 
channels [4] . This technique may have the advantage 
of eliminating the sidelobe issue for certain hardware 
configurations in LAS systems, while allowing for 
multiple channels that share the same bandwidth. 
Additionally, the time -shifting approach allows for 
simple and straightforward data processing of multi- 
channel laser returns since only a single matched- 
filter correlation needs to be performed for all 
channels instead of a separate matched-filter corre- 
lation for each channel. One potential disadvantage 
with this technique is that not all hardware configu- 
rations have a suitable low-frequency response for 
transmitting and/or receiving a pure ML sequence 
without distortion. Since an ML sequence has a main 
power band at zero frequency, using such a technique 
in band-limited hardware (high-pass or bandpass fil- 
tered in particular) will filter out the main power 
band and corrupt the autocorrelation properties, 
which can result in nonorthogonality between chan- 
nels in a time shifted multichannel configuration. 
This problem could be addressed by imparting a 
PN code amplitude modulation on a narrowband 
carrier at a predetermined frequency, performing 
amplitude demodulation, and conducting matched 
filter correlation with the reference PN code to obtain 
the range capability. A low-cost demonstration of 
this was implemented using Hilbert transforms for 
amplitude demodulation. This demonstrated much 
better autocorrelation properties than a linear swept 
frequency but still exhibited minor artifacts, and the 
demodulated PN code exhibited minor distortions 
and ringing [5]. This was a student demonstration 
implemented in audio and was instrumental in 
demonstrating remote sensing concepts using PN 
codes on a carrier and was also a nice demonstra- 
tion outside of the classroom as well. In this paper, 
we eliminate those minor issues by implementing 
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an improved demodulator. We also introduce im- 
proved modulation techniques beyond amplitude 
modulation, which improves the signal-to-noise ra- 
tion (SNR). To further simplify our swept frequency 
IM-CW LAS system for a future ASCENDS mission, 
we concurrently investigate advanced IM techniques 
and develop modulation/demodulation techniques 
that are simple and have autocorrelation properties 
that rival those of the pure ML sequence. This 
includes a quadrature correlator demodulator, PSK 
modulation techniques, and hybrid-pulse modula- 
tion. This results in near-perfect autocorrelation 
properties and an increase in SNR by up to a factor 
of 2. 

Different modulation techniques are discussed in 
the following sections appropriate for the hardware 
shown in Fig. 2. Key results are summarized in 
the conclusion. Details of the actual instrument, 
atmospheric models, etc. have been discussed else- 
where [11,18] and will not be repeated here. 

2. Measurement Technique 

The concept of our IM-CW system for CO 2 IPDA mea- 
surements is illustrated in Fig. 1. For the system de- 
scribed in this figure, amplitude modulated on- and 
offline seed lasers with good spectral properties are 
optically combined using fiber optics and used to 
simultaneously seed a single erbium-doped fiber 
amplifier (EDFA) to increase the transmitted power. 
A small fraction of the transmitted beam is picked off 
via an optical tap inside of the EDFA and sent to a 
reference detector for energy normalization. The 
backscattered science signal from the on- and offline 
wavelengths from the surface as well as aerosols and 
clouds are simultaneously collected with a telescope, 
optically filtered with a narrowband optical filter, 
and detected by a single detector. The science and 
reference signals are amplified, electronically fil- 
tered, and then digitized for retrievals of column 
C0 2 . Postprocessing of the digitized science and 
reference data allow for discrimination between 
ground and intermediate scatterers (from clouds, 
etc.) using the matched-filter technique and also 
to obtain differential absorption power ratios for 


inference of C0 2 column amounts as well as range 
estimates to the scattering targets. 

The on- and offline laser modulation signals using 
the time-shifted approach for this study can be 
expressed as 

A 0 fr = 1 + m£(t),A on = 1 + m£(t - At), (1) 

where m is the modulation index with a value be- 
tween 0 and 1, and £(t) is the repeating modulation 
waveform where -1 < f(f) < 1, {c(t)) — 0, and A t is 
the time shifting of online modulation waveform 
related to that of offline. In general, we choose At > 
2r m ax/c to avoid interference between the time- 
shifted channel and intermediate scatterers, where 
r max is the maximum distance among potential scat- 
terers, including the Earth’s surface, and is also 
called the unambiguous range, while c is the speed 
of light. 

Given the modulation waveforms presented in 
Eq. (1), and using the relation R^ n/off (D = 
■^n/offAon/offfr), the instantaneous received lidar 
on- and offline received optical powers [P^ ff (t) and 
Po n (t), respectively] from a single scattering target 
(i.e., the surface) at range r (which may also vary 
with time depending on the experiment) are given by 

= f ^ exp (-2c /V)dr') 

x exp(-2r) (1 + m%(t -2 r/c)) 

Ponit) =^PLexp(-2£^ /l(r')dr' j 

xexp(-2r)exp(-2r')(l + m%{t-2r/c- At)), (2) 

where t is the total one-way column optical depth for 
the on- and offline measurement, which is deter- 
mined by the gas absorption from molecules other 
than C0 2 , and t is the one-way column optical depth 
resulting from C0 2 absorption only, [’> is the backscat- 
ter coefficient (km -1 sr -1 ), e is the extinction to back- 
scatter ratio (sr -1 ), A is a constant, where P^ ff and 


On-line modulation 



Fig. 1. Baseline instrument block diagram with either bandpass or low-pass-filtered receiver subsystem. Current baseline uses a 
bandpass filter. 
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Pon are the average transmitted on- and offline 
powers averaged over one period of the modulation 
waveform, respectively. 

When multiple targets including the surface exist 
in the lidar path length, each target would generate 
similar received lidar signals like those in Eq. (2). 
The signals from these multiple targets are com- 
bined at the detector and converted to an electronic 
signal S(t). For an AC coupled receiver subsystem, 
we have 

S(t) = Y^ c ik*n£(t - 2 r k /c - At) + C 2k m£(t - 2 r k /c)\, 

k 

(3) 


where we sum over all k scatterers and where, 


c ik= I ^tPln ex v(~ 2T f 0 Vfr'jdr' j exp(-2^)exp(-2^) 

C 2k =^^off ex P (~ 2>: j () V(>')dr'j exp(-2r fe ), (4) 

where K' is a constant. These Cy t and C 2 * returns can 
be uniquely discriminated from other returns from 
different scatterers using a matched filter, which 
involves correlating the received waveform with a 
reference signal. These are also proportional to the 
average received optical power and P[j n k ) from 
the Mh scatterer. For a ground target, solving for z g 
gives 



\c lg p T J 



w’ 


(5) 


where C lg and C 2g are the ground returns. Once C lg 
and C 2g are determined, the column optical depth for 
C0 2 can be found using Eq. (5). In general, this is 
done by cross correlating the reference waveform 
with the return signal S(t), i.e., by the matched-filter 
technique as mentioned before. This results in multi- 
ple peaks for the multiple scatterers including the 
surface. The ground return can be identified through 
range gating (or maximum range). By choosing the 
reference modulation waveforms with perfect auto- 
correlation properties (ML sequence, for instance), 
channel separation is achieved through the time- 
shifting approach since each channel correlates at 
a different apparent range without interference be- 
tween channels or scatterers. That is, each scatterer 
results in two peaks in its matched-filter output: one 
for the offline returns; the other for online returns 
with the At shifted time. The details of this approach 
are discussed elsewhere [4]. 


3. Bandpass-Filtered Modulation Waveform 
Generation and Demodulation 

We now present the improved modulation approach 
for C0 2 IPDA measurements. This new method mod- 
ulates an IM sine wave carrier by an ML sequence as 


in Eq. (6). To effectively represent each bit of ML 
sequence in the transmitted waveforms and allow 
flexibility in the unambiguous range for a particular 
code length and sample rate, we oversample the ML 
sequence by an integer, M, such that each code bit of 
the ML sequence is represented by M points. We then 
use that bit to modulate a sine wave that has exactly 
P cycles per code bit. In order to satisfy the Nyquist 
rule, we must have M > 2 P. Let z n = z(n) be the origi- 
nal ML sequence and Z n = Z{n ) be the oversampled 
ML sequence such that Z(n ) = z(int((n - 1 )/M) + 1), 
where int(x) is a function that represents the integer 
part of x. These definitions are used in the following 
sections wherever PN code modulation is discussed. 
The following two subsections describe amplitude 
and PSK IM carrier modulation and demodulation 
techniques. Amplitude modulation is achieved by 
modulating the IM carrier with an ML sequence. 
With PSK modulation, we use the ML sequence to 
shift the phase of the IM carrier by 180 degrees, 
which is equivalent to changing the sign of the IM 
carrier. The PSK modulation technique makes more 
efficient use of available optical amplifier power and 
results in better SNR. 

A. IM Carrier Amplitude Modulation 

In the case of amplitude modulation, the modulation 

waveform takes on one of two forms: 

£ a (n) = Z(n) cos(2 n(n - 1/2 )P/M), 

^b(n) — Z(n) sin(2^(n - 1/2 )P/M). (6) 

The Vi factor is introduced to make the autocorrela- 
tion function symmetric. 

Demodulation is accomplished as in Fig. 2 using 
the reference kernel of 2 Z{n) - 1. This particular 
reference kernel is necessary to obtain a resulting 
autocorrelation function with off main lobe values 
of zero with the amplitude modulation described 
by Eq. (6). That is to say, with this reference kernel 
the autocorrelation function is zero except at the 
main lobe. This is a quadrature demodulator [ 19 ] 
with the low-pass filter replaced by a matched-filter 
correlation with the oversampled ML sequence. The 


Demodulator 



Fig. 2. Quadrature-matched filter correlator used for demodula- 
tion is similar to a standard quadrature demodulator except the 
low-pass filter section has been replaced by a matched-filter cor- 
relation with the reference PN code kernel. 
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correlations can be calculated using a DFT, which are 
most efficiently computed with a fast Fourier trans- 
form (FFT) and take the form [4,5], 

^ N - 1 

f?(ref , data) = — ^ ref*(m)data(m + n) 

m = 0 

= DFT _1 (DFT*(ref)DFT(data)), (7) 

where * denotes the complex conjugate, N is the 
length of oversampled ML sequence, ref is defined 
as the reference signal, and data is defined as the 
received data. 

As an example of a modulation technique that is 
suitable for a bandpass-filtered LAS system pre- 
sented in Fig. 1, we take the special case where P — 
1 (one cycle per code bit) and M — 8 (8 samples per 
code bit). We generate an 8th order (or 255 length) 
ML sequence using the recurrence relation z n+8 — 
z n ® z n+ 2 ® z n+ 5 ® z n +6 with the seed {1,0, 1,0,1, 
1, 1, 1}, where © is the “exclusive or.” We then over- 
sample this by a factor of 8. Figure 3 shows a com- 
parison between the cosine wave versus sine wave 
modulation as in Eq. (6). Note that each waveform 
in this and all subsequent plots are plotted with eight 
samples per code bit with linear interpolation 
between points. 

The autocorrelation function may be computed 
according to the block diagram in Fig. 2. The autocor- 
relation function for the cosine and sine wave cases 
are shown in Fig. 4, using a 35 sample delay. For 
larger values of P, the autocorrelation for the sine 
and cosine modulation cases approaches that of 
the pure PN code autocorrelation function. 

It can be shown that the half-height width for the 
cosine and sine wave cases is identical to the pure 
ML-sequence case, though the sine wave case may 
be better suited for absorption measurements be- 
cause of the more rounded (flatter) peak. The maxi- 
mum value for an autocorrelation function with a 
sharp peak will show sensitivities to under-sampling 
and therefore may result in a bias. The reason for 
this is because the return signal will be delayed 
in such a way that the exact peak location of the 


autocorrelation function may not always fall on a 
sample point. As a result, there will be a range- 
dependent jitter in the height measurement. If the 
peak is sharp, the jitter will be worse than if it is rel- 
atively flat. Another reason is that typically one uses 
some sort of interpolation to get a better estimate of 
the peak height, and this is much easier to do if the 
peak is a continuous function. The sine wave autocor- 
relation function is less sharp than the pure ML 
sequence. In both cases, and within the accuracy of 
numerical calculations, the off main lobe values 
are essentially zero with the chosen parameters 
and normalization. The resolution of discriminating 
different targets is close to the pure ML-sequence 
case of c/2B r [4,5], where B, is the code bitrate (num- 
ber of code bits per second), which is related to the 
effective modulation bandwidth. It should be noted 
the autocorrelation function in the case of more than 
two cycles per code bit results in an autocorrelation 
function that looks very similar to the pure ML 
sequence case [4,5] . 

B. PSK IM Carrier Modulation 

Our particular PSK modulation takes one of two 
forms: 

C a (n) — (2 Z(n) - 1) cos(2;r(n - 1/2 )P/M), 

? b (n) = (2 Z(n) - 1) sin(2;r(n - 1/2 )P/M). (8) 

As in the previous case, the V 2 factor is introduced to 
make the autocorrelation function symmetric. This is 
very similar to the amplitude modulation case except 
that Z(n) has been replaced by 2Z(n) - 1 Figure 5 
shows a comparison between the cosine wave versus 
sine wave modulation as in Eq. (8). 

Demodulation is accomplished as in Fig. 2, which 
results in the best autocorrelation properties. 
Another PSK demodulation technique sometimes 
used in radar is to correlate the outgoing PSK modu- 
lated sine wave signal with the reference PSK modu- 
lated sine wave, but that results in less than ideal 
autocorrelation properties that contain side lobes 
[20]. In this case, we use the reference kernel Z(n) 
with aforementioned demodulator in Fig. 2. This 



(a) 


— pn code 


(b) 


Fig. 3. Comparison between cosine (a) and sine wave modulation (b) for amplitude-modulation case as in Eq. (6). Each modulation wave- 
form produces a different autocorrelation function. 
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Fig. 4. Comparison between autocorrelation function for cosine wave (a) and sine wave case (b) with ML-sequence modulation with one 
cycle per code bit and eight samples per code bit using a 35 sample delay. Half-height width is identical to the pure ML-sequence case. 


reference kernel is necessary to get a resulting 
autocorrelation function with off main lobe values 
of zero in the case of PSK. The autocorrelation func- 
tion is identical to that shown in Fig. 4 in the case 
where P — 1 (one cycle per code bit) and M — 8 (eight 
samples per code bit). In spite of identical autocorre- 
lation functions, this case has a major advantage. 
Since we are using Z(n ) as a kernel, this method 
has better noise immunity. The reason is because 
Z(n) is 0 for half of the length, which means we cor- 
relate over half the noise compared to the amplitude 
modulated case. As a result, the SNR should be bet- 
ter by a factor of s/2 since the amplitude for each is 
the same. However, if we were to use 2 Z(ri) - 1 for the 
kernel in the PSK modulation case, the SNR would 
be even better. The noise increases by a factor of s/2, 
but the amplitude increases by a factor of 2. The net 
result is that the SNR is a factor of 2 better than the 
amplitude modulation case. These results are ex- 
plained in detail in Appendix A. The downside is 
there is a slight constant bias in the background 
that is inversely proportional to the bit length of 
the code [4], which makes it undesirable for the 
amplitude measurements necessary for C0 2 retriev- 
als. This should present little or no problem for 
range measurements, however. As a result, the 
Z{n) kernel version used here is the likely candidate 


for amplitude measurements because of the superior 
autocorrelation properties. 

4. Low-Pass Filtered Modulation Waveform Generation 
and Demodulation 

A. Hybrid Pulse/CW Modulation 

In lidar designed to detect pulses such as those with 
photon-counting receivers or other types of DC 
coupled receiver subsystems, other types of modula- 
tion are possible. For an alternative implementation, 
we consider the case where we have DC-coupled elec- 
tronics, which, for the purpose of this paper, are 
achieved by using the low-pass-filtered configuration 
in Fig. 1, combined with an EDFA in the transmitter. 
An EDFA amplifier utilizes an internal pump laser 
with a wavelength typically at 980 nm that continu- 
ally pumps the erbium atoms in the erbium-doped 
fiber to an excited state. An external signal laser 
(which in our case has a wavelength of approximately 
1570 nm) is used to pass photons through the 
erbium-doped fiber. Excited erbium atoms then give 
up some of their energy at the same wavelength and 
phase as the signal laser, which has the effect of am- 
plifying the laser signal. There are two types of 
modulation possible with an EDFA amplifier. The 
first is pump modulation. Modulating the intensity 
of the pump laser produces a change in gain of the 



( a ) — pn code 


Fig. 5. Comparison between cosine (a) and sine wave modulation (b) for PSK modulation case as in Eq. (8). Each modulation waveform 
produces a different autocorrelation function. 
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Fig. 6. EDFA optical amplifier step response (signal modulation). 


signal, which modulates the output of the EDFA. 
With pump modulation there is a low-pass-filter re- 
sponse [ 21 ] . This would be an ideal configuration, but 
it requires an optical amplifier for each channel. 
However, with signal modulation as has been used 
in the past [11], which allows one to use the same 
amplifier for all channels, this system behaves as 
a lead compensator circuit [21], which is similar to 
a high-pass filter. This is accomplished by modulat- 
ing the input laser signal directly. In our case, we 
modulate this signal by varying the intensity. The 
step response of such an amplifier configured for sig- 
nal modulation, which is shown in Fig. 6, does not 
fully follow the input step function. This is due to 
an excited state population depletion, which then 
settles to a steady state. 

Even if the transmitter and receiver electronics 
are DC coupled, the nature of the optical amplifier 
used here conspires to degrade the autocorrelation 
properties due to the distortion of the PN code. By 
modulating the PN code by an offset carrier, one 
can potentially avoid those sorts of distortions be- 
cause the optical amplifier responds better to quick 
changes in signal input than it does to maintaining a 
constant level as is required by a pure PN code. One 
solution is to build this as a pulse train into the 
modulation signal. One could do this with a square 
wave if one has enough bandwidth. A particularly 
simple method would be to modify Eq. (6) to 

£(n) = (1 - cos(2;r(n - 1/2 )P/M))Z(n) - 1 

= 2 sin 2 (n(n - 1 /2 )P/M)Z(n) - 1 . (9) 

An example of this modification is shown in Fig. 7 
with the first 100 points of the modulation waveform. 
For this example, P — 1 and M — 8 are selected. This 
modification has one main advantage over the pre- 
vious technique: for the same amplitude, it uses only 
half the average optical power. The reason is because 


— Modulated sine wave 

— pn code 


”'’0 20 40 60 80 100 

Sample number 

Fig. 7. Plot of hybrid sine wave pulse laser modulation signal, 
using modulation waveform shown in Eq. (9), assuming m = 1. 
Note that the received signal in a DC coupled system will be pro- 
portional to the modulation signal. 
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Fig. 8. Autocorrelation function for hybrid pulse modulation us- 
ing a 35 sample delay. 


in order to transmit the previously described modu- 
lation it must be offset. This means that a logical zero 
in that case is half the peak laser power whereas in 
this case a logical zero is zero power. 

Demodulation of this is particularly simple. 
Instead of demodulating the result using the block 
diagram in Fig. 2, one simply has to correlate the 
return signal with 2 Z(n) - 1. The resulting autocor- 
relation function is shown in Fig. 8. 

As in the previous case, the half-height width is 
identical to the pure ML-sequence modulation. How- 
ever, due to the more rounded peak of the autocorre- 
lation function, the height measurement, which is 
desirable for accurate C0 2 absorption measure- 
ments, should be less sensitive to sampling than with 
the sharp peak of the pure ML-sequence case or the 
ML-sequence modulated sine wave case in the pre- 
vious section. Yet the autocorrelation properties 
are near perfect as those in previous cases with 
the chosen parameters and normalization. 

5. Bandwidth-Limited Codes 

In most LAS implementations it is desirable to limit 
the bandwidth of the PN code to reduce aliasing and 
thereby minimize sampling errors in the measure- 
ment of the amplitude height of the autocorrelation 
function. This is particularly important in a space ap- 
plication because there is an upper limit to how fast 
one can sample. For instance, as of the time of this 
writing, the fastest one can digitize with space quali- 
fied parts, at 16 bit resolution is 10 MHz, which es- 
sentially sets the upper limit for the spatial distance 
between points at about 15 m. Since we need 3 m 
resolution accuracy to get 1 ppm C0 2 concentration 
accuracy, it puts severe limits on what we can do 
without interpolation, and interpolation doesn’t 
work well with an aliased signal. Since we can’t sam- 
ple at any rate we want to avoid aliasing, the only 
option left is to somehow band limit the signal. This 
helps with the absorption measurements and makes 
the peak height less sensitive to small changes in 
range. The question is how to accomplish this 
requirement? Many filters can introduce sidelobes 
to the autocorrelation function, such as a square 


window filter or some other high-order filter. This 
also potentially includes a nonoptimal bandpass fil- 
ter, which may already be present so filtering the PN 
code in the correct way before hand in combination 
with putting the PN code on a carrier can prevent 
the problem. Gaussian filters are sometimes used 
in digital communications to bandwidth limit digital 
data, minimize channel bleed-over, and carry the 
maximum amount of information with the least 
bandwidth [22], A convenient implementation is to 
use a periodic version of a Gaussian filter so that 
we can terminate the filter at the end of a cycle. Since 
a discrete Fourier transform (DFT) contains the fre- 
quency content in the first half of the array and a 
mirror image in the second half, and the PN code 
itself is transmitted as a repeating waveform, a cyclic 
filter convolution is appropriate. 

We may construct a periodic Gaussian in the fol- 
lowing way with 


G(x) — 22 exp((x - nL) 2 /2a 2 ). (10) 

n =- oo 

Here L is the period, a is the standard deviation, x is 
an as yet undefined variable, and L is the period. 
This sum can be done analytically (see Appendix B) 
as 


Y, exp((x - nL) 2 /2a 2 ) 


= V2n 

= sl2n 


1 + 2 22 ( ?" 2 cos (2 nnx/L) 


n = 0 

& 3 {nx/L,q), 


(ID 


where $ 3 is the Jacobi theta function of the third 
kind, and q is given by 


q = exp(-2;z-V/L 2 ). (12) 

Theta functions are widely used in the theory of 
elliptic functions [ 23 ] and other applications 
[24,25]. In order to scale this so that it is 1 at the 
maxima, we have 


W(x,L,q) = 8 3 {nx/L,q)/8 3 (0,q). (13) 

To apply this as a filter in the frequency domain, we 
simply multiply the DFT by the array W(n - l.N.q) 
and take the inverse DFT, where N is the length of 
the array and 1 < n < N. Here we have taken x — n - 
1 with period N since the analytic continuation of a 
DFT has period N if n is the array index. 

To demonstrate the usefulness of this filter, we 
first show how to limit the bandwidth of an ML 
sequence and how the autocorrelation properties 
and resolution are preserved. We then use that 
bandwidth-limited sequence in the two modulation 
examples discussed in the previous sections to show 


the effect of bandwidth limiting and preservation of 
the near-perfect autocorrelation properties for each 
case with only a minor degradation in resolution. 
Because of these characteristics, the discussed modu- 
lation codes could be easily implemented in current 
C0 2 LAS systems and the concept of future space 
instruments. 

A. Bandwidth-Limited ML Sequence 

We start with the consideration of the simple case of 
a pure ML sequence with length 255 and over- 
sampled by a factor of 8 as in the previous section. 
We filter this with the theta function filter with q — 
0.9 and N — 255 x 8 = 2040. A plot of the unfiltered 
versus filtered PN code frequency content is shown 
in Fig. 9. 

If the original PN code modulation were applied to 
a LAS system with a nonoptimal filter, aside from the 
spread in the autocorrelation function, the spectrum 
content could be distorted in a way that could 
introduce possible artifacts such as sidelobes. By 
bandwidth limiting the codes, it is possible for the 
modulation bandwidth to fit well within the band- 
width of the bandpass or low-pass receiver filter such 
that the receiver filter impact is minimal to virtually 
eliminate the sidelobe issue, which has been con- 
firmed by simulations. These concepts point out 
practical ways to implement our PN codes to current 
and future C0 2 LAS instrumentation. 

The waveform shapes of these unfiltered versus fil- 
tered PN codes are shown in Fig. 10. The waveform of 
the filtered code is obtained from the circular convo- 
lution between the cyclic theta function filter and the 
cyclic PN code. 

The results of the correlation between 2 Z(n) - 1 
and the filtered and unfiltered PN codes are shown 
in Fig. 11. In spite of the strong filtering, the main 
lobe width of the autocorrelation function is similar, 
and the peak is more rounded, which is a desirable 
trait for accurate measurements of C0 2 column 
amounts. Also, if the waveform is sampled at higher 
than the Nyquist rate, the peak height of the 
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Frequency (sample number) 

Fig. 9. Comparison of frequency content of the filtered versus 
unfiltered PN code. Filtered PN code bandwidth is well under 
the Nyquist sampling rate demonstrating aliasing is no longer 
an issue. 
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Fig. 10. Comparison between filtered PN code and unfiltered PN 
code. 

autocorrelation function can be determined more 
accurately through interpolation because the con- 
tinuous function it is based on can always be recon- 
structed exactly Shannon’s version of the Nyquist 
theorem states exactly that [26]. Of course, noise 
and other errors would tend to make this process less 
exact. But it does show that bandwidth limiting the 
waveform to below the Nyquist sample rate is desir- 
able if one is interested in absorption measurements. 
Because of the required IM carriers for PN modula- 
tion in the instrumentation as discussed previously, 
we further investigate the characteristics of sine 
wave and hybrid modulations for the filtered PN 
codes to ensure the usefulness of current discussed 
techniques for LAS applications. 

B. Amplitude Modulated IM Carrier Using Theta Function 
Filtered PN Code 

In this case, we use the filtered PN code in the pre- 
vious subsection to modulate a sine wave IM carrier 
as in c h of Eq. (6) with P — 2 (two cycles per code bit) 
and M — 8 (eight samples per code bit). A plot of the 
resulting modulation is shown in Fig. 12. 

The frequency spectra of the modulated sine wave 
are plotted in Fig. 13. This shows that the resulting 
frequency content of the filtered modulation is 
bounded within the Nyquist range with no aliasing. 

The autocorrelation functions of these filtered 
versus unfiltered PN modulations are illustrated 



Fig. 11. Autocorrelation function for filtered versus unfiltered 
PN code (filtered case was correlated with unfiltered reference 
using a 35 sample delay). Off-pulse values are zero to within 
numerical precision. 



Fig. 12. Filtered PN code and resulting amplitude modulated 
sine wave. 



Fig. 13. Frequency spectra of the filtered modulation versus un- 
filtered modulation. Modulation bandwidth for the filtered case is 
under the Nyquist sampling rate demonstrating aliasing is not an 
issue. 

in Fig. 14. This shows the modulation produces more 
robust matched filter output for peak determination 
by producing a more rounded peak, which is desir- 
able for amplitude measurements in differential 
absorption estimation because it produces less sam- 
pling error. The resolution degradation is very minor. 



Fig. 14. Comparison of autocorrelation functions between fil- 
tered modulation case and unfiltered modulation case using a 
35 sample delay. Off-pulse values are zero to within numerical 
precision. 
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Fig. 15. Filtered PN code and resulting PSK modulated sine 
wave. 


C. PSK Modulated IM Carrier Using Theta Function 
Filtered PN Code 

In this case, we use the filtered PN code described 
earlier to modulate a sine wave IM carrier as in % b 
of Eq. (8) with P — 2 (two cycles per code bit) and 
M — 8 (eight samples per code bit). A plot of the 
resulting modulation is shown in Fig. 15. 

Figure Ifi is a plot of the frequency spectra of the 
modulated sine wave. This shows the resulting fre- 
quency content of the filtered modulation is also 
bounded within the Nyquist range with no aliasing 
as in the amplitude modulation case. 

The plot of the autocorrelation function of the 
filtered versus unfiltered PN PSK modulation is 
identical to Fig. 14 of the previous section. 

D. Hybrid Pulse Modulation with Theta Function Filtered 
ML Sequence 

We take the case using Eq. (9) with P — 1 (one cycle 
per code bit) and M — 8 (eight samples per code bit) 
and use the same filtered PN code as in the previous 
two examples with q — 0.9. Figure 17 shows the 
resulting modulation. 

The resulting frequency spectra are shown in 
Fig. 18. Here we see that the original modulation 
with the chosen parameters already does a good 



Fig. 16. Frequency spectra of the filtered modulation versus un- 
filtered modulation. Modulation bandwidth for the filtered case is 
under the Nyquist sampling rate demonstrating aliasing is not an 
issue. 



Fig. 17. Hybrid pulse modulation with filtered PN code modula- 
tion. 



Frequency (sample number) 


Fig. 18. Comparison of the frequency spectra between the filtered 
and unfiltered hybrid pulse case. In this case, the filtered PN code 
modulation bandwidth is under the Nyquist rate, but the unfil- 
tered modulation has a limited bandwidth too. 

job at limiting the high frequencies as one might de- 
duce from Fig. 6. This is because there are no sharp 
transitions as there are with the pure ML sequence. 
The filtered version limits the spectra even more. 

Figure 19 is a plot of the autocorrelation function 
of the filtered versus unfiltered PN code modulation. 
Both cases show very good autocorrelation prop- 
erties. 

6. Discussion and Summary 

Table 1 shows a summary of the various modulation 
techniques covered in this study. We have presented 



Fig. 19. Comparison between the filtered versus unfiltered auto- 
correlation functions for the hybrid pulse case, using a 35 sample 
delay. Off-pulse values are zero to within numerical precision. 
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several different PN modulation techniques that can 
be used with different types of hardware: the base- 
line system in Fig. 1 with a signal-modulated EDFA, 
the baseline with the low -pass filter configuration, or 
the baseline system with the low-pass configuration 
and using separate EDFAs with pump modulation. 
In the case of DC coupled electronics combined with 
a laser amplifier with pump modulation, a pure or 
filtered ML sequence not on an IM carrier is ideal. 
In the case of the baseline, although the ampli- 
tude-modulated sine wave version works well, the 
PSK modulation case discussed here has a superior 
SNR. If the electronics is DC coupled, but the laser 
has a transient response such as the case with a sig- 
nal modulated EDFA, the hybrid pulse version is 
best. In addition, a new type of cyclic filter based 
on the Jacobi theta function was introduced, which 
is useful for bandwidth-limiting digital PN codes 
without seriously degrading the resolution and re- 
sulting in autocorrelation properties completely 
without side lobes or other artifacts. The main ben- 
efit is to reduce aliasing, which reduces sampling 
error when measuring the amplitude of the autocor- 
relation function. It can also minimize the impact of a 
nonoptimal filter by minimizing side lobes caused by 
that filter as has been confirmed by numerical sim- 
ulations with filters used in previous experiments. 
By using a combination of IM carrier modulation 
and theta function filtering, one can design a modu- 
lation scheme that is flexible enough to work with 
most hardware while maintaining the integrity of 
the autocorrelation properties. Internal studies have 
shown the previous linear-swept -frequency tech- 
nique, which exhibits a sidelobe level of up to 
13 dB, used in our previous test range and flight cam- 
paigns, has shown a degradation in accuracy due to 
sidelobe interference from thin clouds. Although one 
could potentially develop numerical techniques to 
deal with this issue, it is much simpler to eliminate 
the problem by using a modulation and demodula- 
tion that doesn’t have sidelobes. This helps meet 


Table 1. Comparison of Several Different Modulations* 



Total BW 

PBCF 

SNR 

SB Offset 

PN 

Wide 

DC 

2s/2 

No 

FPN 

Variable 

DC 

<2s/2 

No 

AM 

Wide 

Carrier 

1 

No 

FAM 

Variable 

Carrier 

<1 

No 

PSK1 

Wide 

Carrier 

s/2 

No 

FPSK1 

Variable 

Carrier 

<s/2 

No 

PSK2 

Wide 

Carrier 

2 

Yes 

FPSK2 

Variable 

Carrier 

<2 

Yes 

H pulse 

Less wide 

DC 

2s/2 

No 

FH pulse 

Variable 

DC 

<2sf2 

No 


“Total bandwidth (total BW), power band center frequency 
(PBCF), signal-to-noise (SNR) using AM modulation as a 
reference, and signal background offset (SB offset), 
modulations include pure PN code (PN), filtered PN code 
(FPN), amplitude modulation (AM), PSK with Z kernel 
(PSK1), filtered PSK with Z kernel (PSK1), PSK with 2Z-1 
kernel (PSK2), and filtered PSK with 2Z-1 kernel (PSK2). 


ASCENDS program requirements by reducing com- 
putational burden. Another benefit that has been 
discussed in a previous publication [4] is the use of 
the time-shifted approach to separate channels. That 
approach works with every modulation technique 
presented here and has a major advantage in that 
the demodulation process can be done with a single 
demodulation step for all channels. This further 
simplifies the computational burden when compared 
to the linear-swept-frequency technique, which re- 
quires a separate demodulation for each channel. 
Numerical simulations based on a previous model 
for the system used up to this point [ 18 ] show that 
these techniques could potentially benefit the AS- 
CENDS program with little or no change to the 
existing baseline hardware and will be explored in 
a future experimental study. 

Appendix A: SNR for AM, PSK, PN, and Hybrid Pulse 
Modulation 

Let’s suppose we have a sine wave signal amplitude 
modulated by a PN code, Z with amplitude a in the 
presence of white noise (for example) r] with variance 
a 2 and 0 mean. For simplicity, we look at the pure 
electronic signal as pure signal plus white noise. 
Note that in reality the noise variance will also be 
a function of the average received optical power if 
we have detector shot noise, which is normally mod- 
eled as white noise. If we assume equal average 
power, this simple white noise with variance a 2 
assumption is valid. This is sufficient for determin- 
ing whether amplitude or PSK modulation has better 
SNR, and what it is under those assumptions. 

1. Noise through the Quadrature Correlator 
A. 2Z-1 Kernel Case 

We first look at the noise level at the output of the 
lower and upper correlators in Fig. 2. The variance 
of the correlation of PN code kernel 2Z - 1 with 
the lower sine reference signal multiplied by the 
noise (as per the demodulation process in the lower 
leg in Fig. 2) is 

u|=^((2Z-l) 2 sin 2 (2^> 2 ) 

= ^ (sin 2 (2 nftj)ri 2 ) « ( A1 ) 

so the standard deviation is 

oj k , — - a. (A2) 

L VM 

For the noise at the output of the upper correlator, we 
have 
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(A10) 


4 = i((2Z-l) 2 cos 2 (2^> 2 ) 
= ^(cos 2 (2 nftj)n 2 )~^\o 2 . 
Here the standard deviation is 


Vw a ' 


(A3) 


(A4) 


B. Z Kernel Case 

As before, we look at the noise level at the output of 
the lower and upper correlators in Fig. 2. The vari- 
ance of the noise at the variance of the correlation 
of PN code kernel Z with the lower sine reference sig- 
nal multiplied by the noise (as per the demodulation 
process in Fig. 2) is 

4=^ <Z 2 sin 2 (2 7i ft j)tj 2 ) * ~ol (A5) 


so that 


2 VN 0 ' 


For the upper leg, we have 

4 = ^(Z 2 cos 2 (2^ 7 4)«1^ 2 , 


(A6) 


(A7) 


so that 


“ ^7 N ”■ <A8) 

Generally speaking, the noise at the output of the 
lower and upper correlators {r] L and i]jj) will have 
the same variance but will generally be statistically 
independent because they are both the result of two 
independent correlations. 

2. SNR at the Output of the Quadrature Correlator 

For each modulation, we look at the effect of signal 
and noise through the quadrature correlator. We spe- 
cifically compare the signal and noise at the correla- 
tion peak, since this is how we measure absorption. 
The output of the quadrature correlator is 


Stot — \J (&l + 1 l ) 2 + (Su + flu) 2 ■ (A9) 

It’s clear that the noise level will be reduced by a 
significant amount depending on the length of the 
correlation. We assume the correlation is long 
enough that the correlation peak of the signal is 
much bigger than the noise even if we start with a 
weak signal. Let 


s-^/si + si. 

Then 

Stot = S^i + 2 | \n' L + 2^4 + (|) 2 + (|) 2 
~ 

— S + sin(<j 6 )/ 7 L + cosf^f^j/- (All) 

The output noise is then 

'/out « sin((p)r] L + cosicpiriu . (A12) 

The variance of the output noise is then 

4it = (7out> « ((sin(0)>7L + COS {<p)tiu) 2 ) 

= sin 2 ((j)){rj 2 ) +2 sin( 0 )cos {(p){tl L nu) 

+ cos 2 (^)(4). (A13) 

Since rj L and rju are statistically independent, 
{'Il’Iu) ~ 0- From our previous results, (r/ 2 ) = (r/ 2 ,)- 
Therefore 

4.t = (sin 2 ((p) + cos 2 ((/))) {>il ) =4 = 4 . (A14) 

The signal to noise is 

q q q 

SNR = = — = — . (A15) 

^out n v 


A. Amplitude Modulation Case with 2Z-1 Kernel 

In this case, we use 2Z-1 as the kernel in the quad- 
rature correlator and assume a signal with ampli- 
tude a. The signal peak at the lower and upper 
correlators will be 


S L — sin(</>), S n = cos(^), (A16) 

where cp is the phase angle between the incoming car- 
rier and the sine reference signal. Using the results 
of Eqs. ( A16 ) and (A2), we have 


SNR 


S _ \I S l + s u _ a s/N 

0L °L 2 s[2a ' 


(A17) 


B. PSK Modulation Case with Z Kernel 

In this case, we use Z as the kernel in the quadrature 
correlator and assume a signal with amplitude “a.” 
The signal peak at the lower and upper correlators 
will be 
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S L = -a sin {(p), 

Sir — cos(</>), (A18) 

where cj) is the phase angle between the incoming car- 
rier and the sine reference signal. Using the results 
of Eqs. ( A18 ) and (A6), we have 


SNR 


S _ \I S l+ s u _aVN 

a L o L 2(j 


(A19) 


This is better than the amplitude modulated case by 
a factor of %/2. 


The signal to noise is then 


SNR = — 

^out 


aVN 

a 


(A24) 


This is better than the amplitude-modulation case by 
a factor of 2 \[2. However, unlike the other cases there 
is a bias offset in the noise. The bias of the correlated 
noise is 


rj mt = ((2Z-l)rj)=^fj, (A25) 


C. PSK Modulation Case with 2Z-1 Kernel 

In this case, we use Z as the kernel in the quadrature 
correlator and assume a signal with amplitude “a.” 
The signal peak at the lower and upper correlators 
will be 

S L — sin(0),Sj7 = cos(cp), (A20) 

where </j is the phase angle between the incoming car- 
rier and the sine reference signal. Using the results 
of Eqs. ( A20 ) and (A2), we have 



This is better than the amplitude modulated case by 
a factor of 2. 

D. SNR for the Pure PN Code and Hybrid Pulse Case 

In order to make it comparable sine wave signal case 
for the same average optical power, we take the out- 
put signal of 2aZ. For the case of the hybrid pulse we 
use 4a sin 2 (w(n - l/2)P/M)Z(n) from the results of 
Eq. (9). The extra factor of 2 comes from equalizing 
the average optical power for each (assuming the 
laser amplifier has the headroom to make this pos- 
sible). Using 2Z-1 as the kernel the correlation peak 
for either modulation will be 

N+ 1 

S — — a « a. (A22) 

N 


where N c is the bit length (uninterpolated code 
length) of the PN code. On the other hand, this is 
small compared to the standard deviation. There is 
also a small bias term for the other carrier-based 
modulation cases even if the original noise isn’t 
biased. The small terms we threw away in Eq. ( All ) 
represent a small bias since those noise terms are 
squared. This is a consequence of the quadrature 
demodulation, which makes the net signal plus noise 
always positive. One possible solution to that would 
be to demodulate the signal with a phase-corrected 
reference carrier and use only one leg of the demodu- 
lator. The downside is this may not be simple to do in 
flight conditions where the range and phase are con- 
stantly changing. 

Appendix B: Periodic Gaussian as a Jacobi Theta 
Function 

The sum represented by Eq. ( 10 ) may be performed 
analytically using the Poisson summation formula 
[27], which states 


£/■(«)=£/■(&), (Bi) 

n =- oo k=-oo 

where 


f (k) — / f(n) exTp(-2mkn)dn, (B2) 

J -OO 


The noise in this case will not have 0 mean because 
transmitting a pure PN code usually requires a DC 
coupled receiver or photon counting system. Because 
the detector will output one-sided shot noise, there 
will always be a bias offset in the noise. The variance 
of the correlated noise is 

= I ((2Z - 1) 2 (,, - Tj) 2 ) = ^ (0/ - r}) 2 ) = ^o 2 . 

(A23) 


is the Fourier transform of f. In the case of Eq. (10), 
this becomes 


fik) 


-f 


exp((x - nL) 2 /2a 2 ) exp(-27ri£n)dn 


= s/2n 


exTp(-27t 2 (j 2 k 2 /L 2 )exp(-2nkx/L), (B3) 


so that 
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F(x) — s/2x 


k—-oo 

x exp(-2 nkx/L) 


exp(-2 n 2 a 2 k 2 /L 2 ) 


= V2n 

= \[2n 


[1 + 2 ^2 Q nl cos(2 nkx/L)] 


k=0 

9 3 (nx/L,q). 


(B4) 
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